Proposed experiment: out-of-plane nodal lines in Sr2RuC>4 
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Since the original proposal of an unconventional chiral order parameter in the ruthenate perovskite supercon- 
ductor Sr2RuC>4, much attention has been given to the possibility of out-of-plane nodal lines on the predominant 
y cylindrical Fermi surface given evidence for low-lying quasiparticle excitations in this material. Here I propose 
a tunneling spectroscopy experiment to determine whether such nodal lines in fact exist. 

PACS numbers: 



Introduction. Superconductivity at approximately 1 K was 
discovered in S^RuCU in 1994 by Maeno et al 1 , and since that 
time has been a topic of strong interest, with substantial exper- 
imental and theoretical activity continuing fifteen years after 
its discovery. This material was found following a lengthy 
search for high-temperature materials structurally similar to 
the high-T c cuprates, but not containing Cu. It was thought 2 
that a new family of high-temperature superconductors might 
be discovered in this way, and while this has not happened, 
interest in this material remains high. 

Almost immediately after its discovery, Rice and Sigrist 3 
proposed that this material contained a two-dimensional p- 
wave order parameter as an 'electronic analogue' to superfluid 
Helium. A closely related compound, SrRuC>3, shows ferro- 
magnetism, and so the argument was made in analogy to the 
ferromagnetically mediated pairing in He. 

Complicating this simple picture, however, is substantial 
evidence for nodal excitations in this material. The proposed 
chiral order parameter would give rise to low-temperature 
exponentially activated behavior in the various thermody- 
namic quantities (such as magnetic penetration depth and nu- 
clear spin relaxation rate), but this is not what has been ob- 
served. Bonalde et al 4 measured the London penetration depth 
in single crystals of Si"2RuC»4 and found a T 2 dependence 
,while Nishizaki et al 5 found T 2 specific heat behavior, evi- 
dence for a line-node state. Power-law behavior was also ob- 
served in nuclear spin relaxation rate (T^ 1 ) measurements 6 
and ultrasonic attenuation 7 . In addition, Izawa et al 8 mea- 
sured the magnetothermal conductivity of single crystals of 
Sr2RuC<4 and found that any nodal lines could not be par- 
allel to the c-axis, which immediately suggested A(k) = 
exp(/(j)) cos(c£ z ), given the previous evidence for nodal ex- 
citations. The lack of anisotropy in ab-plane magnetother- 
mal conductivity measurements 9 also suggests any nodal lines 
are parallel to the basal plane. Most recently, Ishida — et al 
again conducted T^ 1 measurements on a high-quality sample 
of Sr2RuC»4 and found T 3 behavior, commonly taken as in- 
dicative of line nodes. Given these measurements, there may 
well be nodes parallel to the basal plane on Sr2RuC»4. 

In this paper I propose an experiment that could help de- 
termine whether these nodal lines exist. The method is based 
on tunneling spectroscopy, which can be a strong probe of or- 
der parameter symmetry. The basis of the experiment is pre- 
sented in Figure 1 (reprinted from 11 ), which depicts, in mo- 
mentum space, an Sr2Ru04 c-axis tunneling spectroscopy ex- 



periment, based upon a recent proposal by the author and P. 
Thalmeier for the use of graphite as a normal electrode in a 
superconducting-insulator-graphite tunneling experiment. As 
the method is described in detail in that publication I only 
sketch the proposal here. The basic idea is that the use of a gate 
voltage applied to the semi-metal graphite changes the length 
of the electron-occupied graphite HKH Fermi surface "cigar". 
When the graphite is deployed in an appropriate c-axis orien- 
tation, as indicated, the conservation of the momentum paral- 
lel to the interface k\\ means that different cigar lengths will 
sample different regions of the Sr2Ru04 Fermi surface. If the 
superconducting order parameter has no k z -dependence, the 
sole effect of the lengthening of the cigar will be an an in- 
crease in conduction channels and thereby merely an overall 
scale factor in the conductance. However, if the order parame- 
ter has k z dependence, each point on the cigar will see a region 
of different A(k), producing tunneling or Andreev density of 
states features at that energy. The differential conductance, 
when properly normalized to the high-bias value, would thus 
vary with gate voltage. 

In this paper I limit myself to a proposal for experimental 
detection of an exp(/0)cos(fc z ) order parameter (OP). Several 
other gap functions have been proposed, including the d-wave 
order parameter exp(/(j))sin(A: z ) 12 , as well as various other d- 
wave and f-wave OP's. The method described herein for the 
exp(/(j)) cos(k z ) will in general yield distinguishable results for 
any OP with significant k, dependence, but due to space con- 
straints I present explicit results only for this OP. 
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FIG. 1 : A proposed tunneling spectroscopy experiment to help deter- 
mine Si"2Ru04 order parameter symmetry. Depending on the gate- 
voltage controlled length of the graphite Fermi HK line, regions 
of different order parameter value are selected, leading to different 
dl/dV behavior. 

Calculation. The calculation presented in this section 
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follows the standard technique applicable to Andreev and 
tunneling spectroscopy of anisotropic superconducting order 



parameters 
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and therefore I keep only the most essen- 



tial details of the calculation. In general, the pair state of 
a superconductor is described by the Bogoliubov-deGennes 
equations 
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with f representing electron-like wavefunctions and g repre- 
senting hole-like wavefunctions, with solutions 



/(x,k,f) = «(k)exp(i (k-r — Et)/h) 
g(x,k,t) = v(k)exp(i(k-r+Ef)/7i) 

with u and v are the BCS coherence factors 1319 : 
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Here (]) is the phase of the gap A(k). Given an electron in- 
cident from the normal metal, two additional particles result 
in the metal: an Andreev-reflected hole 20 , and a normally re- 
flected electron, while in the superconductor an electron-like 
and hole-like quasiparticle result. Later in this work we allow 
for the effect of quasiparticle scattering by letting the energy 
E have a finite imaginary part r 16 i 21 i 22 i 23 . 

Each of the particles above has a corresponding amplitude 
(a, b,c and d, respectively) which is found by specifying the 
boundary conditions: continuity of the wavefunction across 
the boundary, and the following condition applicable to 8- 
function barrier potentials, as introduced ini^: 



V|/ s (0)- ViV (0) = -jHv(0) 
n 
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with the barrier function potential H8(jc). Once a and b, 
the amplitudes for Andreev and normal reflection, have been 
solved for the differential conductance dl/dV is calculated 
from the following: 



dl/dV • 



J dk\\{\ 



\a(k h E)\ 2 -\b(k h E)\ 2 ) 



(8) 



The calculated results employ the boundary condition that k|| , 
the momentum parallel to the interface, is conserved. This 
condition allows for wavevector selection along the longitudi- 
nal Sr2RuC>4 Fermi surface, and hence the acquisition of in- 
formation about the order parameter value at this wavevector. 

We note that for such an ab-plane Si^RuCU tunneling ex- 
periment, the phase of the order parameter is of crucial im- 
portance, as it determines the phase of the hole component 
wavefunction v. Depicted in Figure 2 is an ab-plane diagram 
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FIG. 2: The ab-plane layout of the proposed experiment in 
momentum-space. Note that the real-space graphite axes are rotated 
30° with respect to the momentum space axes. 



of the proposed experiment. The basic point is that assuming 
parallel momentum conservation, the electron-like quasiparti- 
cle will contact the Sr2Ru04 Fermi surface at k/, = 0, k fl = k? 
sampling the order parameter at this wavevector, while the 
backscattered hole will sample the wavevector directly oppo- 
site to this, which introduces a relative phase of -1 among the 
hole-like components of these quasiparticles, due to the order 
parameter sign change. As shown by Tanaka 15 , such a sign 
change leads to a zero-bias conductance peak. However, we 
will see that the character of this peak and surrounding fea- 
tures will depend substantially upon the gate voltage if the 
order parameter has k z dependence, leaving a characteristic 
signature of such a state. 

Main Result. Depicted in Figure 3 is the main result of 
this paper. The differential conductance dl/dV is shown in 
two limits, for several gate-induced graphite chemical poten- 
tial changes for the barrier parameter Z: top, Z=0.3, corre- 
sponding to the point contact regime, and bottom, Z=2, corre- 
sponding to the tunneling regime. 

Substantial effects of the gating are apparent in both plots. 
In the top plot, the low-energy Andreev signal evolves from a 
rounded hump containing a sharp feature at V = Ao towards a 
nearly linear behavior for 8/j = — \2.1meV, while in the bot- 
tom plot the depression in dl/dV that may sometimes occur 
adjacent to a ZBCP, most prominent for 8/j = —20meV grad- 
ually fills in as the gate voltage increases and lower energy 
states are accessed, additionally narrowing the peak itself. 

The evolution of the curves with gate voltage is easily un- 
derstood. For the Andreev-limit plot, the width of the Andreev 
reflection (AR) signal narrows with increasing gate voltage 
because one is seeing an AR signal from a portion of Fermi 
surface with smaller A(k) than the maximum gap, and the 
sharp feature present at V = Ao for 8/j=-20meV similarly be- 
comes less prominent. Analogously, for the tunneling limit 
the width of the ZBCP decreases with increasing gate voltage, 
and the feature at V = Ao is washed out by the summing of 
ZBCP curves with progressively smaller effective A(k). 

These differences are sufficient that a point contact or tun- 
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FIG. 3: The results of dl/dV calculations in the Andreev limit 
(top, Z=0.3) and tunneling limit (bottom, Z=2), for graphite gated 
are shown. Graphite chemical potential changes 8jj = —20meV 
(solid,black), — 18.8raeV (dashed,red), — ll.3meV (dot-dashed, 
blue), — l5.3meV (double-dot-dashed, green) and — 12.7 meV 
(double-dash-dotted, maroon). Insets: Dynes T taken as 0.75 An. 



neling experiment performed along these lines should be able 
to distinguish them. If there is no cos(k z ) dependence to the 
order parameter, the dl/dV curves for these varying gate volt- 
ages should be essentially identical (up to a scale factor). 

For simplicity, we have chosen above the real-space orien- 
tation orientation of the graphite such that the hexagonal face 
parallels the interface. Substantially different dl/dV results 
obtain if instead the hexagonal face is perpendicular to the in- 
terface, as shown in the k-space Figure 2. In this case the 30" 
rotation of the Brillouin zone relative to the real-space unit 
cell means that the graphite Fermi lines no longer occur at 
k[, = 0, but are displaced above and below by approximately 
0.737/A, which is relatively near the Sr2Ru04 y band kp of 
0.75/A, so that the perfect order parameter sign change de- 
scribed above does not occur and one does not see a ZBCP. 
Results for this case are presented in Figure 4, for the same 
parameters as in Fig. 3, and as in the previous plot substan- 
tial gate-voltage created differences are apparent, with the An- 



dreev limit peak at V — Ao reducing with increasing gate volt- 
age, while in the tunneling limit substantial sub-gap density- 
of-states appear with increasing gate voltage. 
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FIG. 4: The results of dl/dV calculations in the Andreev limit (left, 
Z=0.3) and tunneling limit (right, Z=2), for graphite gated and ro- 
tated 30" relative to that of Fig. 3 are shown. Same parameters as in 
Fig. 3. 

It is clear that these methods remain applicable even if the 
k z dependence of the order parameter is a + bcos(k z ), with 
\a\ < b. The dl/dV curves will still evolve with gate voltage. 
The main effect of such a constant term would be to change 
the gate voltage at which the Fermi surface cigar accesses the 
nodal lines, with a minor secondary effect on the shape of the 
DOS, so long as a <§; b. 

Discussion - Experimental Consideration: Surface degra- 
dation. The results of the above section suggest that appro- 
priate Andreev or tunneling experiments may be able to de- 
termine whether or not there are nodes in the gap function 
located at k z = ±Jt/2, as an exp (/(])) cos (£,) order parameter 
would contain. There is, however, an experimental considera- 
tion requiring consideration: the surface degradation effects 
in this material. 

There are to date three spectroscopic measurements made 
on Sr2Ru0 4 25 ' 26 . The first, by R. Jin et al 26 , performed tun- 
neling on a Pb-Sr2RuC>4 junction but found only a spectrum 
strongly resembling that of superconducting lead, with coher- 
ence peaks (assumed to be that of Pb) at approximately 1.4 
meV and no sub-gap structure. The second experiment, by 
Upward et al22 performed c-axis STM using a Pt/Ir tip and 
found clear evidence of a superconducting gap. The third ex- 
periment, by Laube et al 27 , observed a weak Andreev reflec- 
tion signal, as well as a zero-bias anomaly. Of issue for our 
proposed experiment is the extremely large zero-bias conduc- 
tance - approximately 0.85 the normal-state value in the Up- 
ward et al data. This is believed to result from surface degra- 
dation. 

A detailed accounting for the possible effects of surface 
degradation on tunneling or Andreev spectra is rather compli- 
cated, particularly as the cause of such degradation remains 
unknown. Rather than attempt a first principles calculation as 
such, we therefore content ourselves with two relatively sim- 
ple mechanisms for simulating the effects of surface degrada- 
tion. 
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FIG. 5: dl/dV curves assuming significant nonconservation of paral- 
lel momentum, as described in the text. Same parameters as figure 3, 
top (left here) and bottom (right). 

The first method is empirically based upon the comparative 
smallness (i.e. relative to background conductance) of the fea- 
tures observed in the spectroscopic work performed to date, 
and the lack of detail in these features. Both of these char- 
acteristics are consistent with a large effective "smearing" of 
the conductance. To a lesser degree, such smearing is nearly 
universally present in spectroscopy on unconventional super- 
conductors, where it is typically modeled by adding an imag- 
inary part r to the quasiparticle energy 21 . In fact such a T 
has already been taken to represent the effects of a disordered 
surface layer 28 . Given the predominance of this "smearing" 
in experimental spectra of unconventional superconductors, 
and the comparative success of the Dynes T in modeling such 
smearing, it is reasonable to use this parameter to simulate the 
as-yet-unknown source of surface disorder in S^RuOzt. 

I have performed a calculation, for the graphite orientation 
used in the main panels of Fig. 3, with a large Dynes T of 
0.75Ao; results of this calculation are presented in the insets 
of Fig. 3. Even with this large T, the figures continue to 
show substantial differences in the dl/dV curves. The zero- 
bias conductance changes significantly from 8/j = —20meV 
to — HJmeV in both the Andreev and tunneling regime, as 
does the shape and size of the Andreev reflection signal, so 
that even if the problems with surface degradation persist, the 
nodal lines, if existent, will leave their distinctive signature in 



dl/dV. 

A second method of simulating the effects of surface disor- 
der is to relax the assumption of perfect parallel momentum 
conservation. The premise here is that the surface degradation 
causes scattering within the interface region, so that a finite, 
rather than infinitesimal region of superconductor Fermi sur- 
face is allowed to receive current from a given graphite Fermi 
surface location. In practice one may use a Gaussian distri- 
bution of conductance (i.e. « (1 + \a\ 2 — |Z?| 2 )exp(— a(k\\ — 
fcrjj ) 2 > where k ,|| represents the parallel momentum of the in- 
coming electron and k|| represents the parallel momentum of 
the transmitted quasiparticles, and observe how the dl/dV fea- 
tures change. We note that the boundary conditions of the 
Calculation section do not explicitly involve the parallel mo- 
mentum, so this scheme is internally consistent. The results 
are shown in Figure 5, where we show several Andreev limit 
dl/dV curves assuming an a of 1/V0.2jt (i.e. Ak\\ ~ 0.2ji) 
in units where 2n is 2k/ c, c the Sr2Ru04 lattice constant of 
12.72 A. Here Sfcii is of the same order as typical kt\ ~ 7t/2, 
so that the non-conservation of parallel momentum is substan- 
tial. We have also allowed for non-conservation of k|| in the 
kft direction, i.e. perpendicular to the graphite Fermi line, and 
taken account of the A(k) phase factors that result. Even in 
this case, significant differences in dl/dV remain, so that the 
proposed experiment can be considered robust against the sur- 
face degradation issue. 

Conclusion. In this work I have demonstrated that a series 
of superconductive NIS tunneling measurements on Sr^RuCU 
using gated graphite as the normal electrode should allow de- 
termination of the presence or absence of order parameter line 
nodes parallel to the ab-plane. I have shown that the power 
of such an experiment to determine nodal structure is robust 
against the surface disorder prevalent in this material. I await 
the results of such experiments with great interest. 
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